Transition from Poissonian to GOE level statistics in a modified Artin's billiard 
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' One wall of Artin's billiard on the Poincare half plane is replaced by a one-parameter (cp) family 

I of nongeodetic walls. A brief description of the classical phase space of this system is given. In 

^ . the quantum domain, the continuous and gradual transition from the Poisson like to GOE level 

' statistics due to the small perturbations breaking the symmetry responsible for the 'arithmetic 

( , chaos' at Cp = 1 is studied. Another GOE Poisson transition due to the mixed phase space for 

. large perturbations is also investigated. A satisfactory description of the intermediate level statistics 

' by the Brody distribution was found in both cases. The study supports the existence of a scaling 

, region around Cp = 1. A finite size scaling relation for the Brody-parameter as a function of 1 — Cp 

• and the number of levels considered can be established. 

0\ ■ 

■ 05.45+b,03.65.-w 

V ■ 

O i I. INTRODUCTION 

O ' One of the fundamental questions of quantum chaos 
' is the relationship between the different types of classi- 
• i-H I cal motion and the energy level statistics of the quantum 
^ ' counterpart of the system ||J|,§1 . One important field for 
^ i testing ideas on the subject is the geodesic flow on con- 
- - I stant negative curvature surfaces y,^. The most often 
studied cases are the symmetric octagon and the modular 
surface [@,§j7|J^ . The level spacing statistics of these sys- 
tems were not found to follow the predictions of the ran- 
dom matrix theory for the Gaussian orthogonal ensem- 
ble (GOE) [Ijep. This is contrary to the general belief 
that completely chaotic systems with time reversal sym- 
metry have GOE level spacing statistics [pl flO . This 
unexpected feature was explained recently pUpl based 
on semiclassical arguments. Unlike generic ergodic sys- 
tems, there is an exponential degeneracy of time periods 
of closed classical orbits caused by nontrivial sym- 
metries [ pdl , which is responsible for the amplification of 
the quantum interference in these systems f^. 

In this paper we investigate a modification of Artin's 
billiard introduced in Ref. j7|, where one geodetic wall 
is replaced with a one parameter family of nongeodetic 
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walls. The advantage is that changing only one parame- 
ter we can study perturbations of this billiard. The per- 
turbation breaks the degeneracy of time periods and the 
transition from Poissonian to GOE level spacing distribu- 
tion can be followed. At parameter values corresponding 
to large perturbations another transition from GOE to 
Poissonian statistics can be observed due to the change 
of the ergodic part of the mixed phase space up to com- 
plete integrability p^ . 

The paper is organized as follows. In Sec. |l| we in- 
troduce the family of the billiards, and give a brief de- 
scription of the classical dynamics. In Sec. |l| we shall 
investigate the manifestation of the classical dynamics in 
the statistical properties of the energy levels. Finally, 
Sec. 1^ is devoted to conclusions. 



II. CLASSICAL TREATMENT 

Our billiard is given in the Poincare half plane with 
Gaussian curvature K = — I and Riemannian scalar cur- 
vature R — —2. The motion between two bounces can 
be obtained from the classical Lagrangian 

L=^(i'+f) (1) 
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or Hamiltonian 



where 



H^ly'{pl+Pl) (2) 



X y 

Px = — , Py^ — • 3 

y V 

The domain of the classical motion is the area bounded 
by the vertical walls 

a; = 0, x=l/2 (4) 

and the circle 



X 



2 ' (y-l + i?)2 = i?2 (5) 



with radius R, and center (0, 1 — R). We use the inverse 
radius Cp = as a control-parameter. The collisions 
with the walls are elastic. The trajectories of the free 
motion between the bounces are circles centered on the 
X axis. The footpoint {x) and the inverse radius (M = 
1/r) of the trajectories (Fig. |l|) is a good choice for the 
"bounce mapping" T. The pair {x,M) remains constant 
during the free motion and it changes into a new pair 
(x', M') at each bounce 

(a;,M) -> (x',M') = T{x,M). (6) 

We have studied the classical dynamics in the parame- 
ter range < Cp < 2. Note that for Cp > 2 the bottom of 
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the billiard is open. In the classical treatment we mea- 
sure the time in (2£^)~^/^ units. Then the elapsed time 
agrees with the metric length of orbits. The billiard at 
Cp = 1 is the domain of Artin's triangle billiard and is 
bounded by three geodetic walls. At Cp = the billiard 
has a horizontal non-geodetic straight segment, and the 
motion is integrable [ p^ , p^ . 

In the parameter range 1 < Cp < 2 the motion is com- 
pletely chaotic, the Lyapunov exponent is positive for 
each periodic orbit as we will show next. For Cp = 1 
this is a consequence of Artin's proof for the ergodicity 
[ psf . For other parameters we have used a method com- 
mon in the mathematical literature for the 
calculation of the Lyapunov exponent of the trajectories. 
Namely, following the time evolution of the local curva- 
ture K{t) of an orthogonal curve along a trajectory and 

choosing k(0) > the quantity exp( K(t)dt) measures 
the stretching orthogonal to the orbit . Consequently, 
the Lyapunov exponent of the orbit is given by the aver- 
age value of K 



A= lim ^ / K{t)dt. (7) 



r 

T^oo T , 

During the free motion K,{t) is described by the Riccati 



equation |17 



kit) = -K{t)- K^{t) (8) 

where K{t) is the local Gaussian curvature, which has 
a constant value of K{t) = i^T = — 1 in our case. The 
change of k before and after the collision with a wall 
having curvature q in the place of the bounce, is given by 

2q , , 

= K_ (9) 

COS ip 

where Lp is the angle of the bounce. If the curvature q is 
negative or zero for each wall, it is easy to see from the 
solution of (^) that after sufficient but finite time {T„iax) 
K becomes larger than one and then K(t) > +1 holds 
for t > Tjnax ■ The average of k also becomes larger than 
which means that the system is hyperbolic. This is 
the case when Cp > 1, since the curvature of the bottom 
wall is 

where Vg is the radius of the geodesic touching the wall 
in the point where the curvature is computed (Fig. 1.). 
The quantity (|ic| ) is negative or zero, if 1 < Cp = 
The Lyapunov exponent is A = 1, if Cp = 1, and A > 1, 
if Cp > 1. Since these arguments hold for each orbit, the 
system cannot contain stable periodic orbits. 

In the remaining parameter range < Cp < 1 it is not 
obvious for which parameter values the system is com- 
pletely chaotic, and where the first island of stability is 
created. The bottom wall is locally focusing in this pa- 
rameter range pl|. We expect that the highly unstable 
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free motion suppresses the focusing on a global scale up 
to a certain strength of the local focusing. In other words, 
there exists a parameter c* < 1 so that for c* < Cp < 2 
the system is completely chaotic. 

To have an insight how the mixed phase space sets in 
we have investigated analytically the stability of a family 
of periodic orbits starting at the left corner of the billiard 
and bouncing m times on the vertical walls (Fig. ||). This 
family of orbits exists in the whole < Cp < 2 parameter 
range. They are obviously unstable for Cp > 1. We shall 
show later that for Cp < 1 their stability criterions give a 
good insight into the structure of the mixed phase space. 
Varying the parameter, the shape of these orbits remains 
unchanged. The solution of the Riccati equation (||) for 
Cp < 1 and < k(0) < 1 between two collisions with the 
bottom wall reads as 

K„(t) =tanh(t-4'"^) (11) 

where tl^^ is a parameter to be adjusted. According to 
(a) bouncing with the bottom wall, the curvature changes 



Km(i + 0) ==K„(t-0)-Q. (12) 

For the family of the orbits considered here Q takes the 
form 

cos Lpm 

where 1 — Cp is the curvature of the bottom wall in the 
bounce point, and ipm is the collision angle of the orbit 
labelled by m. A direct calculation leads to the equation 

to/2 , , 

cos(^„ — ===. (14) 
^1 + m2/4 

Since according to (|^) the Lyapunov exponent requires 
the average of K{t) only in the limit when T goes to infin- 
ity, one can calculate this with the asymptotic solution 
of K{t) which satisfies the periodicity condition 

K„i{t + T,n) ^ Km{t), (15) 

where Tm is the period of the orbit given by 



sh(T,„/2) = ^l + myA (16) 



from simple geometric considerations. Solving (15) with 
the condition (^2|) leads to the equation 

tanh(-4™^) = tanh(r„ - 4"^) - Q (17) 

for f^. This is a second order equation for tanh(to). Its 
relevant solution is 

tanh(4"^) = ^(-y'c2(TO2+4)-4cp 

+ (1 -Cp) v/4 + TO2y (18) 
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The Lyapunov exponent is the average of k for one pe- 
rfod: 



1 I 

X,n = Tf- / tanh(t-4"Vi 

Jo 

= ^ln(cosh(r™ - 4™Vcosh(4"))). (19) 

From Eq. (|l|) and one can calculate the Lyapunov 
exponent as a function of Cp and the index of the orbit 
m analytically. Fig. ^ shows the Lyapunov exponent 
as a function of Cp for this family. For Cp = 1 all the 
orbits have the same Lyapunov exponent (Am = +1) as 
a well known consequence of the geodetic walls. Using 
Eqs. ( p^ ) and (18) one can also solve the equation 



A„,(cp)=0 (20) 

and get that parameter value Cp where the m-th orbit of 
the family gets stabilized when Cp is decreased 



Cp = 7- (21) 

First the shortest orbit (to = 1) becomes stable at Cp — 
0.8. 

To decide wether there exists another periodic orbit, 
which is stable for a larger parameter value Cp, we have 
investigated the bounce mapping (||) (see Fig. ^. We 
have not found any island in the range 0.8 < Cp < 1. The 
existence of a stable orbit causing a very small island, 
of course, cannot be excluded, but the phase space in 
this range can be considered completely chaotic at least 
practically. 

After this description of the classical domain, we turn 
to the quantum treatment of the system. 



III. LEVEL STATISTICS AND SCALING 
PROPERTIES 

The energy eigenvalues of the billiard can be computed 
by solving the Schrodinger equation 

1 2 fd^^J d'^ilj 



-2^l^ + #j-^^ ^''^ 

with Dirichlet ('0 = 0) boundary condition on the walls 
(^ and (H). In the integrable case (cp = 0) the partial 
differential equation (|2^) is separable in x and y, and the 
eigenf unctions can be derived jl^]. Using this functional 
basis for the non-integrable problem the Schrodinger 
equation (^2|) can be solved numerically. For further de- 
tails we refer to our earlier papers on this model |l|,|l3|,§. 
The unfolding of the levels has been performed by the 
transformation 

x^ = N{Ei), (23) 
where N{E) is given by the Weyl formula of Rcf. 0. 
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Once one has the unfolded eigenenergies, one can ana- 
lyze their fluctuations. We concentrated on the distribu- 
tion p{s) of spacings between nearest neighbour unfolded 
levels (s = Xi+i—Xi), since this has been found to be very 
sensitive the classical dynamics [p2PJlO|]. The distribu- 
tion p{s) and its first momentum are normalized due to 
the unfolding (p^): 

/•oo poo 

/ p{s)ds = 1, / sp{s)ds = 1. (24) 
Jq Jo 

It is expected, that the low lying energy levels are de- 
termined by the quantum nature of the system while 
the higher semiclassical part of the spectrum reflects the 
classical properties, and exhibits universal behavior for 
integrable and completely chaotic systems. Generic inte- 
grable systems show Poissonian distribution p^ , p^ , p5| , p6[ 

p{s) = e-^ (25) 

In the completely chaotic case the spectra can be de- 
scribed by ensembles of random matrices, namely the 
spacing distribution of the GOE applies This dis- 

tribution is almost identical with the Wigner surmise 

Pi^) = f (26) 

In the case of mixed phase space the situation is less 
clear, the fluctuations are non-universal, and the discus- 
sion about relevant parameters characterizing the level 
spacing distribution of these systems is not finished yet. 

It is of particular interest that at Cp = 1 the level 
statistics do not follow the generic behavior of classically 
chaotic systems. In this 'arithmetic' case there 
exist arguments for the relevance of the Poissonian dis- 
tribution based on the Selberg trace formula and on the 
asymptotic distribution of time periods of periodic or- 
bits [^. However, these considerations are valid for the 
asymptotic, infinitely high part of the spectrum. In the 
statistics of the first say 2000 levels deviations from the 
Poissonian statistics can be observed |@,||. For parame- 
ters slightly different from Cp = 1 the universal GOE dis- 
tribution should show up immediately in the asymptotic, 
infinitely high part of the spectrum, if the universality of 
GOE for completely chaotic systems is true. The change 
from Poissonian to GOE will gradually extend towards 
the lower lying levels, when further decreasing Cp. 

According to the classical behaviour of the system and 
the considerations above, we can expect three regions in 
the parameter space where the nearest neighbour level 
spacing distribution of a finite number of levels behaves 
differently. The first region is < Cp < 0.8, where with 
the exception of the integrable case Cp = the phase 
space is mixed one can expect some intermediate statis- 
tics between GOE and Poissonian. The second region is 
0.8 < Cp < 2, Cp 7^ 1 ± A, where GOE is expected. 
The narrow region around Cp — 1, indicated by A is a 
third region, where we should observe certain deviations 
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from the GOE for any finite number of levels due to the 
'arithmetic chaos' in Cp = 1. 

The first and third regions require a familiy of distribu- 
tions interpolating between the Poissonian and the GOE 
for quite different reasons. There are several attempts to 
find a theoretically well established transition formula for 
systems with mixed phase space ||2^,^ . One of the main 
differences in these cross-over formulas is in the small s 
behaviour of p(s). The Berry- Robnik formula gives 
a constant value for small s (p(s) const), while the 
Lenz-Haake formula j29j leads to a linear level repulsion 
s {p{s) ~ s). 

There are other formulas which interpolate between 
( p6| ) and ( psf ) on a pure empirical base. One of them 
was proposed by Israilev and another one by Brody 
Since they are not necessarily associated with a 
mixed phase space they can be tried also in the third 
region. Both of them show power law level repulsion 
for small s {p{s) ~ s'^). To characterize intermediate 
distributions we have chosen the one parameter family of 
Brody distributions since it is simpler than the other one 
and, for other purposes, we could not see any advantage 
of Israilev's distribution. Moreover, in a recent study 132] 
both of them were found equally satisfactory. The one 
parameter family of Brody distributions is given by 

p{s, (3) = as^^ exp(-6s'3+i) , (27) 



where a and b can be derived from the conditions (|2j) as 
a function of the Brody parameter (3 



fi + 2 



a=(/3 + l)6, 6= (r( ^ ) ) . (28) 



For /3 = (Eq) and for /3 = 1 (26) are recovered. The 



cumulative level density of the Brody distribution is 

I{s)= [ p{x)dx ^l-eyip{-bs>^+^). (29) 
Jo 

We have fitted Brody distributions to our empirical I{s) 
functions applying the method of least squares. Using 
the first n energy levels at a certain Cp and finding the 
best fit yields an empirical (3n{cp) value. Note that we 
do not claim, that the Brody distribution is the correct 
asymptotic form, but only that the fitted /3 parameter 
is a good measure where the intermediate distribution is 
lying between ( psf ) and (p6|). We have carried out this 
program in the parameter region < Cp < f.2, calcu- 
lating about 2000 energy levels. We have found that 
the Brody distrribution provides a satisfactory fit in the 
whole parameter region 

In Fig. ^, where the Brody parameter /3 as a function 
of Cp is depicted, one can see that (3 increases continu- 
ously with Cp in the region < Cp < 0.8, reflecting the 
fact that in the phase space the volume of the chaotic re- 
gion grows up. This is in complete agreement with earlier 
findings of Ref. in another system with mixed phase 
space. The saturation of P starting at Cp ~ 0.5 indicates. 
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that the level spacing distribution is not sensitive to small 
stable islands. A relatively sharp transition from GOE 
toward Poissonian can be observed in the 0.1 >| 1 — Cp | 
parameter range, which we have studied in more detail. 

In the critical point Cp = 1 the numerical level spacing 
distribution is still far from the theoretical Poissonian dis- 
tribution. To see the convergence of the Brody parameter 
towards (3{cp = 1) = we have plotted /3„(cp = 1) as a 
function of the first n energy levels (Fig ^.). The conver- 
gence is quite apparent and /?„ can be well approximated 
in the range n = 200 — 2000 as 

/3„(cp = 1) = 13.05n-". (30) 

with a ~ 0.55. This is very close to a ~ ^/V^ depen- 
dence, which is typical in random matrix description of 
transitions between ensembles . When Cp = 1 — S < 1 
we observed for sufhciently small S values {5 < 0.01) the 
following behavior of /? as function of n. For small n 
numbers {n < n*{S)) the /3n(cp = 1 — 5) curves coincide 
and they tend to constant values for n > n*{6) within 
the n < 2000 range (Fig. These values are scaling 
with 6 as 

/3„(1 -6)^6^ (31) 

where 7 w 0.44. The separate scaling relations ( ^0[) and 
( ^ ) can be unified into one finite size scaling relation 
familiar from the theory of phase transitions 

.^—-=S-y{l/5n% (32) 

Pn(l - 0) 

where this relation defines the scaling function /(x). For 
X — > the scaling function goes to a finite value, while 
for X ^ c» it behaves like f{x) ^ x~'^ . The n depen- 
dence ( |30| ) can be recovered by choosing e — a/7. The 
existence of such a formula is a stronger evidence of the 
scaling then the two scaling relations ( |30| ) and ( pl| ) to- 
gether, since the scaling functions computed from differ- 
ent data must coincide. This has been checked at three 
different S parameters. Satisfactory agreement with the 
scaling hypothesis (|3^) has been found within the statis- 
tical fluctuations of the scaling functions (see Fig. ^). 

Since our billiard system is expected to be completely 
chaotic for small 6 values a finite (3^1 asymptotic value 
would mean that the GOE statistics is not universal. We 
think instead, that the finite value of (3 is restricted to 
a region n*{6) ^ n <C n**{S) only. The emergence of 
such a region can be understood in terms of semiclassical 
considerations. The strength of the perturbation of the 
arithmetic billiard is proportional to 6 = 1 — Cp. Let 
us denote by Sp{E) the action belonging to the periodic 
orbit p. The perturbation destroys the degeneracy of the 
actions Sp{E). If two actions in the arithmetic billiard 
are degenerate Sp{E) = Spi{E), p ^ p', they differ in an 
amount of 

ASpyiS,E) = Sp{S,E) ~ Sp,{S,E) = i2E)'/^lp,p,{S) 

(33) 
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in the perturbed case, where lp.p'{S) = lp{6) — Ipi (S) is the 
difference of the metric length of the periodic orbits p and 
p' ■ lp.p'{S) is independent of the energy. The absence 
of the degeneracy of the periodic orbits p and p' can 
be observed in physical quantities when the difference of 
actions is comparable to h (in our units h = 1): 

ASp^p>i6,E)^h. (34) 

Because of the energy dependence ( |3^ ) of this relation, 
at low energies (small n) the absence of the degeneracy 
of actions cannot be observed. The condition (^4|) gives 
a critical energy value to each pair of periodic orbits, and 
their minimum defines n* Emin/d, where d is the mean 
level spacing 0. This argument also explains that all 
the /3n(l — S) curves coincide for n ^ n* . For 6 = this 
region extends to the infinite n range. When (5 > for 
sufficiently high energy (large n) the broken degeneracy 
of all periodic orbits can be observed and the system 
looks like an ordinary completely chaotic system. This 
defines another scale n** depending on Cp. Between n* 
and n** we expect to observe roughly a constant value of 
Pn- From n** , the /3„ function is supposed to converge 
toward = 1. For Cp —^ I n** will grow to infinity, 
and the length of the scaling region is increasing too. 
However, for Cp sufficiently far from 1 the scaling region 
is very restricted or disappears (n* « n**). At parameter 
6 = 0.03 instead of a plateau, /?„ has merely a minimum 
around 700 (Fig.||). 

IV. CONCLUSIONS 

In this paper we have analyzed the connection of clas- 
sical dynamics and the statistical properties of the quan- 
tum energy spectrum in a constant negative curvature 
billiard family with a nongeodetic wall. In cases, where 
the phase space is mixed, we have found remarkable 
agreement with the Brody distribution which interpo- 
lates between GOE and Poissonian. The Brody parame- 
ter was found to be a monotonously increasing function 
of the chaoticity of the system as measured by the magni- 
tude of the positive Lyapunov exponent. Using the Brody 
parameter we could describe a GOE Poissonian tran- 
sition around Cp = 1 due to the arithmetic chaos. We 
observed a finite size scaling relation among the Brody 
parameter, the number of energy levels and the param- 
eter 1 — Cp, and argued that the existence of the scaling 
region can be explained by the energy dependence of the 
condition (|^ . We hope, the detailed quantitative theory 
of this scaling behavior can be worked out in the future. 
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FIG. 1. Our billiard model on the Poincare half plane. 

FIG. 2. Periodic orbits bouncing once in the left corner and 
m times on the vertical walls. 

FIG. 3. Lyapunov exponents as a function of the Cp pa- 
rameter, for the periodic orbits of Fig. 2.. Only nonzero 
Lyapunov exponents are plotted. 

FIG. 4. The phase space of the bounce mapping (^) for 
Cp = 1, 0.5, 0.07, are shown. Notice, that there is no island 
of stability for a certain M{cp) >| A/ | strip. 

FIG. 5. The Brody parameter as a function of the control 
parameter Cp. Each point, denoted by a + sign was calculated 
from the first 1947 quantum levels. 

FIG. 6. The Brody parameter as a function of the number 
of levels on a log-log plot at parameter Cp = 1 corresponding 
to the arithmetic case. The dashed line represents the best 
power-law fit relation /3„ = 13.05n~°'®®. 

FIG. 7. The Brody parameter as a function of the number 
of levels. The curves are taken at Cp = 1.000, 0.996, 0.995 and 
0.990. The bigger the Cp the lower lying the corresponding 
curve. 

FIG. 8. The scaling function f{x) reconstructed from 
l3n{cp) at parameters Cp — 0.990,0.995 and 0.996. Notice 
that the curves coincide within their fiuctuations. 

FIG. 9. The Brody parameter as a function of the number 
of levels at Cp — 0.970. Notice the verry narow scaling region 
around n — 700. 
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